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Abstract — This work aims at combining adaptive protocol 
design, utility maximization and stoctiastic geometry. We focus 
on a spatial adaptation of Alotia wittiin ttie frameworlc of ad tioc 
networlfs. We consider quasi-static networks in which mobiles 
learn the local topology and incorporate this information to adapt 
their medium access probability (MAP) selection to their local 
environment. We consider the cases where nodes cooperate in a 
distributed way to maximize the global throughput or to achieve 
either proportional fair or max-min fair medium access. In the 
proportional fair case, we show that nodes can compute their 
optimal MAPs as solutions to certain fixed point equations. In 
the maximum throughput case, the optimal MAPs are obtained 
through a Gibbs Sampling based algorithm. In the max min 
case, these are obtained as the solution of a convex optimization 
problem. The main performance analysis result of the paper is 
that this type of distributed adaptation can be analyzed using 
stochastic geometry in the proportional fair case. In this case, 
we show that, when the nodes form a homogeneous Poisson point 
process in the Euclidean plane, the distribution of the optimal 
MAP can be obtained from that of a certain shot noise process 
w.r.t. the node Poisson point process and that the mean utility can 
also be derived from this distribution. We discuss the difficulties 
to be faced for analyzing the performance of the other cases 
(maximal throughput and max-min fairness). Numerical results 
illustrate our findings and quantify the gains brought by spatial 
adaptation in such networks. 

I. Introduction 

Stochastic geometry has recently been used for the analysis 
and performance evaluation of wireless (ad hoc as well as 
cellular) networks; in this approach, one models node locations 
as a spatial point process, e.g., homogeneous Poisson point 
processes, and one computes various network statistics, e.g., 
interference, successful transmission probability, coverage (or, 
outage) probability etc. as spatial averages. This often leads to 
tractable performance metrics that are amenable to parametric 
optimization with respect to network parameters (node density, 
protocol parameters, etc.). More precisely, this approach yields 
spatial averages of the performance metrics for given network 
parameters; then the parameters can be chosen to optimize 
performance. This approach takes a macroscopic view of the 
network with the underlying assumption that all nodes in the 
network have identical statistical characteristics. 

In practice, due to randomness and heterogeneity in net- 
works, nodes need to adapt to local spatial and temporal condi- 



tions (e.g., channel conditions and topology) to reach optimum 
network wide performance. For example, nodes in wireless 
LANs adjust their window sizes based on acknowledgment 
feedback; in cellular networks nodes are scheduled based on 
channel conditions and adapt their transmit powers based on 
the measured SINRs, which in turn depend on the transmit 
powers set by other nodes. In all such scenarios, distributed 
adaptive algorithms are used to reach a desired network wide 
operating point e.g. that maximizing some utility. While the 
behavior of such distributed optimization protocols is often 
well understood on a given topology, there are usually no 
analytical characterizations of the statistical properties of the 
optimal state in large random and heterogeneous networks. 

The main aim of this work is to use stochastic geometry to 
study spatial adaptations of medium access control in Aloha 
that aim at optimizing certain utilities. While we identify a 
utility for which stochastic geometry can be used to compute 
the spatial distribution of MAP and the expected utility, we are 
far from being able to do so for all types of utilities within 
the a-fair class and we discuss the difficulties to be faced. 

Let us start with a review of the state of the art on Aloha. 
Wireless spectrum is well known to be a precious and scarce 
shared resource. Medium Access Control (MAC) algorithms 
are employed to coordinate access to the shared wireless 
medium. An efficient MAC protocol should ensure high 
system throughput, and should also distribute the available 
bandwidth fairly among the competing nodes. The simplest 
of the MAC protocols. Aloha and slotted Aloha, with a 
"random access" spirit, were introduced by Abramson [Tl and 
Roberts [17] respectively. In these protocols, only one node 
could successfully transmit at a time. Reference |i4J modeled 
node locations as spatial point processes, and also modeled 
channel fadings, interferences and SINR based reception. This 
allowed for spatial reuse and multiple simultaneous successful 
transmissions depending on SINR levels at the corresponding 
receivers. All the above protocols prescribe identical attempt 
probabilities for all the nodes. Reference [jSJ further proposed 
opportunistic Aloha in which nodes' transmission attempts are 
modulated by their channel conditions. 

Among the initial attempts of MAP adaptation in Aloha, 
reference flT] analyzed protocol model and proposed stochas- 



tic approximation based strategies that were based on re- 
ceiver feedback and were aimed at stabilizing the network. 
References [4], [5\ also optimized nodes' attempt probabili- 
ties (or thresholds) in order to maximize the spatial density of 
successful transmissions. Reference (13] analyzed both plain 
and opportunistic Aloha in a network where all the nodes 
communicate to one access point. They assumed statistically 
identical Rayleigh faded channels with no dependence on 
geometry (i.e., no path loss components). They demonstrated a 
paradoxical behavior where plain Aloha yields better aggregate 
throughput than the opportunistic one. Reference [14] also 
studied optimal random access with SINR based reception. 
However, they considered constant channel gains. They de- 
veloped a centralized algorithm that maximizes the network 
throughput, and also an algorithm that leads to max-min fair 
operation. Reference ifTsl modeled network as an undirected 
graph and studied Aloha under the protocol model. They 
designed distributed algorithms that are either proportional 
fair or max-min fair Reference llTZl built upon the model 
of H, and formulated the channel access problem as a non- 
cooperative game among users. They considered through- 
put and delay as performance metrics and proposed pricing 
schemes that induce socially optimum behavior at equilibrium. 
However, they set time average quantities (e.g., throughput, 
delay) as utilities (or costs), and concentrated on symmetric 
Nash equilibria. Consequently, in their analysis, dependence 
on local conditions vanishes. 

In none of the above Aloha models, nodes account for both 
wireless channel randomness and local topology for making 
their random access decisions, as we do in the present paperj^ 

There is a vast literature on the modeling of CSMA by 
stochastic geometry which will not be reviewed in detail here. 
The very nature of this MAC protocol is adaptive as each node 
senses the network and acts in order to ensure that certain 
exclusion rules are satisfied, namely that neighboring nodes 
do not access the channel simultaneously. However, CSMA as 
such is designed to guarantee a reasonable scheduling, not to 
optimize any utility of the throughput. The closest reference to 
our work is probably [6] where the authors study an adaptation 
of the exclusion range and of the transmit power of a CSMA 
node to the location of the closest interferer This adaptation 
aims at maximizing the mean number of nodes transmitting 
per unit time and space (while respecting the above exclusion 
rules). This mean number is however only a surrogate of the 
rate. In addition, the adaptation is only w.r.t. to the location 
of the nearest interferer. 

We study spatial adaptation of Aloha in ad hoc networks. 
The network setting is described in Section HI] We consider 
quasi-static networks in which mobiles learn the topology, and 
incorporate this information in their medium access probabil- 
ity (MAP) selection. We consider the cases where nodes are 
benevolent and cooperate in a distributed way to maximize the 
global network throughput or to reach either a proportional fair 

'in view of this distinction, we refer to the spatial Aloha protocol of 0) 
as plain Aloha. 



or max-min fair sharing of the network resources. We analyze 
the case where nodes account only for their closest interferers, 
for all nodes in a given ball around them, or even all nodes 
in the network. 

Section HIl] is focused on the distributed algorithms that 
maximize the aggregate throughput or lead to max-min fair- 
ness in such networks. In the proportional fair case, we show 
that nodes can compute the optimal MAPs as solutions to 
certain fixed point equations. In the maximum throughput case, 
the optimal MAPs are obtained through a Gibbs Sampling 
based algorithm. In the max min case, the optimal MAPs are 
obtained as the solution of a convex optimization problem. 

Section |IV] contains the stochastic geometry results. The 
model features nodes forming a realization of a homogeneous 
Poisson point process in the Euclidean plane. We compute 
the MAP distribution in such a network in the proportional 
fair case using shot noise field theory. To the best of our 
knowledge, this distribution is the first example of successful 
combination of stochastic geometry and adaptive protocol 
design aimed at optimizing ceratin utility function within 
this Aloha setting. We also show that the mean value of 
the logarithm of the throughput obtained by a typical node 
can be derived from this distribution. Finally, we discuss the 
difficulties to be faced in order to extend the result to other 
types of utilities. 

The numerical results are gathered in Section |V] The aim of 
this section is two-fold: 1) check the analytical results against 
simulation and 2) quantify the gains brought by adaption 
within this setting. 

II. Network Model 
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Fig. 1 . A snapshot of bipolar MANET with Aloha as the medium access 
protocol. The diamonds represent transmitters, and the connected circles the 
con'esponding receivers. The solid diamonds represent the nodes that are 
transmitting in a slot. 



We model the ad-hoc wireless network as a set of trans- 
mitters and their corresponding receivers, all located in the 
Euclidean plane. This is often referred to as "bipole model" |j3] 



Chapter 16]. There are A*' transmitter-receiver pairs commu- 
nicating over a shared channel. The transmitters follow the 
slotted version of the Aloha medium access control (MAC) 
protocol (see Figure [T]i. A transmitter, in each transmission 
attempt, sends one packet which occupies one slot. Each trans- 
mitter uses unit transmission power We assume that each node 
has an infinite backlog of packets to transmit to its receiver. 
The Euclidean distance between transmitter j and receiver i is 
rji, and the path-loss exponent is a (a > 2). We also assume 
Rayleigh faded channels with hji being the random fading 
between transmitter j and receiver i. Moreover, we assume that 
the random variables h.ij,l < i < N,l < j < N are indepen- 
dent and identically distributed with mean l//x0Thus all hjiS 
have cumulative distribution function (cdf) F(x) = 1 — e^''^ 
with a; > 0. All the receivers are also subjected to white 
Gaussian thermal noise with variance w, which is also constant 
across slots. We assume that a receiver successfully receives 
the packet of the corresponding transmitter if the received 
SINR exceeds a threshold T. 

Let ci be the indicator variable indicating whether trans- 
mitter i transmits in a slot, and pi be is medium access 
probability. Thus P(ei = 1) = pi. When node i transmits, the 
received SINR at the corresponding receiver is 



li = 



Then the probability of successful reception can be calcu- 
lated as follows. 



'{'y^>T\{{hJ,,e,)■.j^i}) 



■ exp ( -mT W^C'jh^t 



Thus 
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where ^ji = t ( ) • Further simplifying. 
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n 



Pj 



1 + 6 



3^ 



(1) 



^The independence assumption is justified if the distance between two 
receivers is larger tlian tlie colierence distance of tlie wireless channel [3]. 
We assume this to be the case. 



Then, the rate or throughput of transmitter i is given by piQi. 

The thermal noise appears merely as a constant multiplica- 
tive factor in the expression for the successful transmission 
probability (see ([T]i). Moreover, in interference limited net- 
works, the impact of thermal noise is negligible as compared 
to interference. We focus on such networks, and thus we ignore 
the thermal noise factor throughout. 

III. Adaptive Spatial Aloha and Fairness 

In this section, we analyze adaptations of spatial Aloha 
that maximize aggregate throughput or achieve proportional 
fairness or max-min fairness. 

A. Maximum Throughput Medium Access 

The throughput maximizing medium access probabilities 
solve the following optimization problem. 



maximize O 



subject to 



<Pi < 1, ieAf. 



l + b 



3 1 



We first argue that the optimum in the above optimization 
problem is attained at one of the vertices of the hypercube 
formed by the constraint set. To see this, suppose p* G [0, 1]^ 
is an optimal solution, and p* e (0, 1) for some i G N. Clearly, 



dpi lp=p* 
= 0. 



1 + fej. 



n 



Pk 



1 + 



Since the partial derivative is independent of pi, pi can be set 
to either or 1 without reducing the value of the objective 
function. This proves our claim. In the following we focus 
only on such extreme solutions. Then the above problem is 
equivalent to finding an M C N such that pi = 1 if and only 
if i G yVf is an optimal solution. Thus we are interested in 



maximize 



E n 



1 + bj. 



An iterative solution: We can pose this problem as a 
strategic form game with the users as players [16 1. For each 
player its action lies in {0, 1}, and the utility function : 
a R is given by 



Uj(0,a_i) = 0, 
Uj(l,a_i) = 1 



n 



jeM\{i} 



1 - 



l + b. 



E 



n 



1 + bit 

jeM\{i} keM\{ij} 



l + b, 



This is a potential game with the above objective function as 
the potential function |15|. Thus the best response dynamics 
converges to a Nash equilibrium. This algorithm can be im- 
plemented in a distributed fashion if each node i knows bij , bji 
for all j, and also M and YlkeM\{j}i^ ^ (1 + bkj)~^) for all 



j & M after each iteration. However, a Nash equilibrium can 
be a suboptimal solution to the above optimization problem. 
To alleviate this problem, we propose a Gibbs sampler based 
distributed algorithm, wherein each node i chooses action 1 
with probability 



Ui(l,a_i)/T 



Pi = 



1 _|_ g«i(l,a„i)/T ■ 

The parameter r is called the temperature. The Gibbs sampler 
dynamics converges to a steady state which is the Gibbs 
distribution associated with the aggregate throughput and the 
temperature r |9|. In other words, we are led to the following 
distribution on the action profiles: 

7rr(a) = w^'eA/''.(a)^ 

where u is a normalizing constant. When r goes to in an 
appropriate way (i.e., as l/log(l + t), where t is the time), the 
distribution nri-) converges to a dirac mass at the action profile 
a* with maximum aggregate utility if it is unique. Notice 
that the aggregate utility J2ieJV^ii^) equals the aggregate 
throughput. Thus the action profile a* is a solution to the 
original throughput optimization problem. 

Remark 3.1: The first two terms in the utility function 
Ui(l,a_i) can be seen as "selfish" part of user i, whereas 
the last summation term is "altruistic" part. The user makes 
a decision based on whether the "selfish" part dominates or 
viceversa. 

Remark 3.2: In a quasi-static network where topology con- 
tinuously changes, although at a slower time scale, different 
sets of nodes are likely to be scheduled to transmit under 
different topologies. Thus, in terms of long term performance, 
maximum throughput medium access is not grossly unfair 



Observe that user i's optimal attempt probability is indepen- 
dent of the attempt probabilities of other users. In particular, 
if fi{l) < 1, user i can perform iterations p'^^^ = fi{Pi) 
autonomously. Furthermore, 



+ bi. 



Pi 



Clearly, \ fi{pi)\ < 1, i.e., fi{-) is a contraction. Thus the fixed 
point iterations converge to the optimal pi starting from any 
G [0,1]. 

Remark 3.3: The characterization of the optimal attempt 
probabilities reflects the altruistic behavior of users. More pre- 
cisely, user i's attempt probability is a function of {bij, j / i} 
which are measures of i's interference to all other users. 
In particular, if 57" < i-^-. if *'s transmission does 

not cause significant interference to the other users, then i 
transmits in all the slots. Unlike the throughput maximization 
problem, there is no "selfish" component in the decision 
making rule. 

Remark 3.4: As the target SINR T ^ 00, 6^- for all 
and the proportional fair attempt probabilities satisfy 

iV- 1 

- = > z = 

P: 



p, ^ 1 - 



Pi 



Pi 



for all i e Af. This yields Pi = for all i £ Af. This 
is expected, because in the limiting case a transmission can 
succeed if and only if there is no other concurrent transmis- 
sion. This is hence Aloha without spatial reuse, and it is well 
known that in this case, the optimal access probability is 1 /N 
asymptotically [7J. 



B. Proportional Fair Medium Access 

The proportional fair medium access problem can be for- 
mulated as follows. 

maximize log {pi qi ) , 

i 

subject to < Pi < 1, i G N . 
The objective function can be rewritten as 

We thus have a convex separable optimization problem. The 
partial derivative of the objective function with respect to pi 

is 1 1 

^ " 5 1 + b„- ~ ' 



Oij - Pi 



which is continuous and decreasing in pi over [0, 1]. We 
conclude that at optimality, for each 1 < i < TV, 

^^^\hip^)---[T.,^rTTrk^y if/.(i)<i, 

1 otherwise. 



C. Max-min Fair Medium Access 

Our analysis in this section follows ifTsl . ifTol . The max-min 
fair medium access problem can be formulated as 



maximize 



subject to < Pi Y\ ( ^ ^ Y 



+ h. 



where constraint functions are defined for all p G [0, 1]^. 
The following is an equivalent convex optimization prob- 
lem (see [18 1 for details): 

1.2 



minimize 



2 ' 



subject to e< \ogp^ + ^ log h - YTb" ) ' * ^ 
The Lagrange function of this problem is given by ["HI 



2 



+ ^Aje-iogp,-^iog(i-^)), 



with Aj > 0, i G A/" being the Lagrange multipliers. 

Minimization of the Lagrange function (which is concave 



in p and 6) gives 



Pi 




if iE,^.^>i, 



otherwise. 



(3) 
(4) 



Wang and Kar I fTSl suggest that the Lagrange multipHers be 
updated using the gradient projection method. More precisely, 
for all i e N, 



A»(n + 1) 



A,(n) + /?(n) 



(5) 



(6) 



where /3(n) is the step size at the nth iteration. Further 
more, [;18, Theorem 2] implies that a solution arbitrary close 
to an optimal solution can be reached via appropriate choice of 
step sizes. However, all the users need to exchange variables 
in order to perform updates. 

Finally, the directed link graph corresponding to our net- 
work is a directed graph in which each vertex stands for a 
user (i.e., a transmitter-receiver pair) in the network. There is 
an edge from vertex i to vertex j in the directed link graph 
if transmission of user i affects the success of transmission of 
user j. Two vertices i and j are said to be connected if either 
of the following two conditions hold: 

1) there is an edge from i to j or viceversa, 

2) there are vertices = i,vi, . . . ,Vn^i,vn = j such that 
vm and v„i+i are connected for m = 0, 1, . . . , n — 1. 

Clearly, the directed link graph for our network model is a 
complete graph; for any pair of vertices i and j there is an edge 
from i to j and also from j to i. In particular, the directed Unk 
graph is a single strongly connected component |'19l. Thus f 19' 
Corollary 1] implies that the above optimization also obtains 
the lexicographic max-min fair medium access probabilities 
that yield identical rates for all the users. 

D. Closest Interferer Case 

Note that a user needs to know the entire topology, and in 
a few cases, also needs to communicate with all the nodes to 
implement the adaptation rules developed in Sections IIII-Al 
IIII-CI In this section, we carry out analysis under the simpli- 
fying assumption that the aggregate interference at a receiver 
is dominated by the transmission from the closest interferer 
This is a reasonable approximation in a moderately dense 
network, specifically when the path loss attenuations are high. 
Throughout this section, we use the notation 

c(i) := argminrji, 
C{i) := {j : c(j) = i}, 



for all 1 < j < A^. In words, c{i) is the strongest interferer of 
node i, and C{i) is the set of nodes to which node i is the 



strongest interferer We assume that there is always a unique 
c{i) for each i. Then, accounting only for the nearest interferer, 
the approximate probability of successful transmission for 
node i is 

Pcji) 



= 1 



l + b, 



The analysis of Sections |III] can be adapted to this simplified 
scenario. 

1 ) Maximum Throughput Medium Access: The throughput 
maximization problem can now be posed as follows. 

maximize O :— ^ ] piQi, 

i 

subject to Q <Pi< N, i £ U. 

As in Section IIII-AI we can argue that some p* G {0, 1}^ 
attains the optimal throughput. Again, an equivalent optimiza- 
tion problem is 



maximize 

MiZAf 



^ / t{c[{)eM} \ 



or alternatively. 



maximize > 1 — > — ; — - 

-^CA/- '^A l + hij 



We now formulate a strategic form game among users, with 
action sets {0, 1} and utility functions given by 



iii(0,a_i) =0, 
Mj(l,a_i) =1 



c{i) 



°2. 

. 1 + 1 



1 + h„(.i\i ^-^ 1 + fe,-,- 



Again a Gibbs sampler based algorithm yields the optimal set 
of the transmitting users. Also, user i only needs to know the 
distances of user c(i) and all the receivers in C(i) and their 
actions to make its decision. 

Remark 3.5: Notice that user i must choose = 1 if 
1 



1 



1 + 6, 



'c{i)i 



^ 1 + 1 



> 0. 



Such users can set their actions to 1, and need not undergo 
Gibbs sampler updates. 

Discussion: Consider a scenario where a node's closest 
interferer does not transmit, i.e., has zero attempt probability. 
Nonetheless, this node always has an active closest inter- 
ferer (unless there are no other nodes in the network). A better 
approximation of the success probabilities, and hence of the 
throughput, is obtained by always accounting for the closest 
active interferer Towards this, let us define 

c{i, A4) :— argmin rji, 
C{i,M):= {jeM:c{j,M)^i}, 



for all i e M. We are now faced with the following optimiza- 
tion problem. 



maximize - , - , 

Mc^^ I ^ 1 + 6j, 

leM \ jec{i,M) ■' / 

We can now define users' utility functions as follows. 

Ui(0,a_i) =0, 
Ui(l,a_i) =1 



c{i,M)i 



E 



1 + bij 1 + 6cO\A4)j 



]<EC{i,MU{i}) 

where M ^ {j e : j ^ i, aj ^ 1} . The analogous distributed 
algorithm (Gibbs sampler based) can again be shown to lead 
to the optimal solution. 

2) Proportional Fair Medium Access: We now aim to solve 
the following optimization problem. 

maximize log {pi qi ) , 

i 

subject to < Pi < 1, i £ N ■ 
Following the discussion in Section IIII-BI we obtain 



otherwise. 



Again, if Ejec(j) ^ > 1' iterations p*^+i = fi{p^) converge 
to the optimal pi starting from any pi e [0, 1]. 

Remark 3.6: If J^jszcii) F~ ^ i-^-' i^ ^'^ transmission 
does not cause significant interference to the users for whom 
i is closest interferer, then i transmits in all the slots. The 
same user may not transmit (in any slot) under the maximum 
throughput objective if 1 — -^^^^ ^ and c{i) and users in 
C(i) transmit. 

We can have explicit formulae for the attempt probabilities 
in a few special cases. 



1) Suppose C(i) is singleton for each i. If C(i) = {j}, then 



Pi 



2 



if bij < 1, 
otherwise. 



2) Linear Network Topology: We now consider a scenario 
where nodes are placed along a line, say R, and are 
indexed sequentially. We also assume that for any node i 
the potential interferers are the two immediate neighbors 
i — 1 and j + 1. This also amounts to assuming C(i) — 
{i — 1, i + 1} for all i. Then, assuming -j^^—^ + ^. "'"^^ > 1, 
the proportional fair attempt probability of node i satisfies 



+ 



1 + h,. 



1 + - Pi 



This is quadratic equation in pi, which on solving give^ 



2 + + 



(f)...-i-&.,.+i)^ 

+ (l+b..i-l)(l+b..i+l) 



3) Max-min Fair Medium Access: Similar to Section IIII-CI 
the max-min fair medium access problem is 



maximize 9, 



subject to 6* < Pi I 1 — 



Pc{i 



1 + 6, 



Again, the constraint functions are defined for all p G [0, 1] 

Let us recall the definition of the directed link graph 
associated with the network. In this section, we only account 
for the interference due the closest interferer Thus there is an 
edge from vertex i to vertex j if and only if j e C{i). We 
assume that the directed link graph is connected (i.e., all the 
vertices in the graph are connected with each other). If it is not 
connected, the max-min fair medium access problem on the 
entire graph decomposes into separate max-min fair medium 
access problems on each of the connected subgraphs, which 
can be solved independently. 

We now pursue the following convex optimization problem 
which is equivalent to the above max-min fair medium access 
optimization problem. 



mimmize 



IE 



subject to 9.1 < log Pi + log I 1 — 



Pc{i) 



1 + 6, 



€JV, 



<9j,j ec(i),j = 



The last set of constraints along with the connected assump- 
tion (of the directed link graph) forces 9i to be equal for all 
i £ J\f. This confirms equivalence to the initial optimization 
problem. Now the Lagrange function is 



ieM leAf \ 



log Pi - log 1 



Pc{i) 



1 + 6, 



+ E E f^iji^t 

iGN" j£C{i)U{c{i)} 



where A,; > 0,fiij > 0,j € C{i) U {c{i)},i e Af, are the 
Lagrange multipliers. An approach similar to Section IIII-CI 
prescribes the following update rules. For all i £ M,j £ 



The other root is greater than 1, and thus is not a valid probability. 



C{i) U {c{i)}. 



= < 



1 otherwise, 

(7) 

-^i- I]jec(i)u{c(i)}(Mij ~ Mji) 

if + J2jec{i)u{c{i)}it^ij ~ A^ji) > 0' 
otherwise, 

(8) 



A,(n + 1) 



A,H + /?(n) ft 



log P i - log 



1 - 



Pc{i) 



1 + bc{i)i 

fi,^{n + 1) =[M,j(n) + I3{n){e, - 6^)] + , 



(9) 
(10) 



where is the step size at the nth iteration as before. 

Observe that any user i can perform updates (l7li-(fT0ll via local 
information exchange. More precisely, it needs to communi- 
cate only with user c(i) and the users in C(i). 

E. A Note on Distributed Implementation 

The aim of the present paper is primarily of theoretical 
nature and it is beyond our scope to discuss implementation 
issues. Let us however stress that the discussed adaptations of 
the MAP are implementable. We will focus on the proportional 
fair case in view the main focus of the paper 

Assume each receiver has a distinctive pilot signal with fixed 
power Since we assumed a quasi-static network in which the 
nodes move at a slower time scale, each node can then learn 
the distance that separates it from a given receiver by listening 
to its pilot signal and by performing a time average over (so 
as to smooth out fading). Once this data is available for all 
receivers, a given transmitter can then solve the key fixed point 
equation that characterizes its optimal MAP. Distinctive pilot 
signals can be obtained by a collection of orthogonal codes 
chosen at random by the receivers. In practice, it is enough 
for a transmitter to detect the "dominant" receivers (i.e., those 
within a certain distance to it), so that the scheme will work 
in an infinite network with a finite (properly tuned) number of 
such codes. 

IV. Stochastic Geometry Analysis 
A. Network and Communication Model 

We now assume that the transmitting nodes are scattered on 
the Euclidian plane according to a homogeneous Poisson point 
process of intensity A. For each transmitter, its corresponding 
receiver is at distance t-q in a random direction. The traffic 
and channel models are the same as in Section As before 
the transmitters use slotted Aloha to access the channel, and a 
receiver successfully receives the packet from its transmitter if 
the received SINR exceeds a threshold T. Finally, transmitters 
adapt their attempt probabilities as described in Section HII] 

Each transmitter is associated with a multi dimensional 
mark that carries information about the adaptive transmission 



probability and the transmission status. Let l> = {Xn,Zn} 
denote a marked Poisson point, where 

' • $ = {Xn} denotes the Poisson point process of inten- 
sity A, representing the location of transmitters in the 
Euclidean plane. 
• {Zn = {<j>n,Pn, £„ )} denote the marks of the Poisson point 
process l> , which consist of three components: 

- {4>n} denote the angles from transmitters to receivers. 
These angles are i.i.d. and uniform on [0, 2n] and 
independent of $. We will call them the primary marks. 

- {pn} denote the MAPs of the nodes; pn is a secondary 
mark (i.e. a functionals of $ and its primary marks, 
see below). 

- {e„} are indicator functions that take value one if a 
given node decides to transmit in a given time slot, 
and zero otherwise. Clearly, P(en = 1) = p„ = 1 — 
P(e„ = 0). In particular, given p„, e„ is independent 
of everything else including {em}m^n- 

The locations of the receivers will be denoted by = 
{Yn = Xn + iro,(t)n)} with (ro,<;/>) ■- (rocos.;/), rosin <;/>). It 
follows from the displacement theorem |2| that is also 
a homogeneous Poisson point process of intensity A. 

The above assumptions will be referred to as the Poisson 
model. We will also consider below a more general case where 
the above marked point process is just stationary. 

B. Proportional Fair Spatial Aloha 

1) MAP distribution: Let us consider response functions 
L : X r2 R+ defined for each < p < 1 as follows 

P 



Lp{x,y) = 



For all < p < IjX G R^, the shot noise field J$r(p, 
associated with the above response function and the marked 
point process l> is 



J$r(p, a;) 



Notice that this shot noise is not that representing the inter- 
ference at X. It rather measures the effect of the presence of 
a transmitter at x on the whole set of receivers. 

Consider a typical node at the origin, Xq = 0, with marks 
Po,(t>o- Let P° denote the Palm distribution of the stationary 
marked point process $ [2, Chapter 1]. The fixed point 
equation determining the MAP of node Xo = reads (see (|2]l) 

We have a similar equation for each node and the sequence 
{pn} is readily seen to be a sequence of marks of $ and its 
primary marks. 

It follows directly from monotonicity arguments that 



1 1 

Po P 



iff 




< 1 



Notice that we have not used the specific assumptions on the 
point process so far. Hence we have the following general 
connection between the optimal MAP distribution and the shot 
noise J$r: 

Theorem 4.1: For all stationary marked point processes 
<& (not necessarily Poisson), for all < p < 1, 



and 



P"(P0 > P) = P° {j^r\{Yo}M < l) , 
P°(P0 = 1)=P° {j<I>r\iYo}ihO)<^), 



with the Palm distribution of 3>. 

We now use the fact that l> is an independently marked Poisson 
point process |2 Definition 2.1]. From SUvnyak's theorem ||2] 
Theorem 1.13], 

P" (^*-\{Yo}(P=0) < l) =P(J$.(p,0) < 1) 

for all < p < 1. Consequently, 

F°(P0>P)=P(J*'-(P,0)<1), 



and 



'(PO = 1)=P(J*K1,0)<1). 



It follows from ||2] Proposition 2.6] and from the fact that 
is a homogeneous Poisson point process that one can write 
the Laplace transform Cj^p o^{s) of the shot noise J<i,r(p, 0) as 



'C,/(p.o)(s) = exp|-27rAy - e '■°+(i-p)-o J rdr 

'(11) 

where fo := Tvq. 

Theorem 4.2: Under the above Poisson assumptions, the 
attempt probability of the typical node has the distribution 



(po > p) = 



2-K 



'Cj(p,o)(«w)- 



-dw, (12) 



with /:j(p,o)(-) given by (dB- 

Proof: Let gp{-) denote the density of the shot noise field 
J$(p, 0). Then 



P°(PO > P) = / 9p{t)dt = / gp 



where u{t) = lifO<t<l and otherwise. Now using 
Parseval's theorem 

1 f°° 

with = Eexp(— iw^) the Fourier transform of the real 

valued random variable A and B* the complex conjugate of 
B. The claim follows after substituting Tu{w) = ^—j:^ — and 
^J(pfi)iw) = Cj(pfi){iw). ■ 
Remark 4.1 : For a = 4, the Laplace transform Cj(^p Q-^{s) 
can be simplified as 



C 



= exp 



Jo V' 



-spv^ I (1— p) 



-d« 



2) Mean Utility: This subsection is devoted to the analysis 
of the mean value of the logarithm of the throughput of the 
typical node: 

e = E°log((po))+E"log((go)), 

with go defined in ([TJ. Since we know the cdf / of po^ the 
first term poses no problem. The second term can be rewritten 
as (see ([B) 



lE°log((go))=EO 



Pn 



||-y„-yo|| 
TVT" 



+ 1 



Under the law P°, the points {Xnjn^o of $ form a homoge- 
neous Poisson point process of intensity A. However, the marks 
{Pn}n^o do not have the law identified in the last section. In 
fact, the mark p„ of a point Xn (n / 0) satisfies the following 
modified fixed point equation: 

1 1 v-^ 1 



— T?°^ 



1 - Pn 



E ■ 



TfT 



+ 1 — Pn 



with the convention that p„ = 1 if there is no solution in 
[0,1]. We can use the same argument as above to conclude 
that J- < i iff 



||x„-yo||° 



+ 1-P 



+ E ■ 



\\Xr.-Yr„ 



< 1. 



+ 1 



Conditioned on there being two nodes at and x, the other 
points form a homogeneous Poisson point process of intensity 
A. This allows one to prove the following. 

Theorem 4.3: Under the above Poisson assumptions, given 
that there is a node at and a node at a; e R^, the attempt 
probability of the node at x has the distribution 

1 f°° e™ — 1 
P"'''(px > p) = TT ^J.(p.o) dw, 



with 



'Cj,(p^O)(s) = 



|a:-(ro,0)||" + (l-p)ro 



exp{-2TT\ I (1-e '■" + (i~p)*'o I r-dr- 



and (rQ,<j)) :— (rocos(j!>, rosin(j!>). 

Due to the circular symmetry, the first integral in the expres- 
sion of £j^(p o)(s) depends on x only through ||a;||. Thus the 
density of px also depends on ||x|| only, and it will be denoted 
by fr when \\x\\ = r. The density of po identified in the last 
subsection will be denoted by /. The main result of this section 



Theorem 4.4: Under the above Poisson assumptions, in the 



proportional fair case, the mean utility of a typical node is 

e = f\og{u)f{du) 

Jo 

^ 2^ I I /^°^ llx-(ro,9i)ll"+fo) 

d<^/||^_(^^_0)l|(dt;)da;. (13) 
Proof: See Appendix lAl ■ 

C. Discussion of the Other Cases 

A preliminary concern when trying to use Euclidean 
stochastic geometry of the infinite plane in the maximum 
throughput and max-min fairness cases is that it is not clear 
whether the associated infinite dimensional optimization prob- 
lems make sense in the first place. In the proportional fair 
case, each node computes its optimizing MAP in one step as 
the solution of a fixed point equation that is almost surely 
well defined (in terms of a shot noise) even in the infinite 
Poisson population case. Unfortunately, this does not extend 
to the other two cases. 

This does not mean that there is no hope at all. Consider 
for instance the maximum throughput case accounting only 
for the closest interferer, and further simplify it by measuring 
interference at the transmitter rather than at the receiver. In 
other words, consider the same optimization problem as in 
Section IIII-Dll but with c(i) being the closest transmitter 
to transmitter i (rather than to receiver t). Then, in the 
Poisson case, the infinite dimensional optimization problem 
can be shown to reduce to a countable collection of finite 
optimization problems. This follows from the fact that there 
are no infinite "descending chains" in a homogeneous Poisson 
point process ifTol Chapter 2] (a descending chain is a sequence 
of nodes ii, 12, . . . such that c(i„) — in+i for all n). As a result, 
the Poisson point process can be decomposed into a countable 
collection of finite "descending trees", where each path from 
the leaves to the root is a descending chain. The associated 
optimization is hence well defined and the problem can be 
reduced to evaluating the solution of the optimization problem 
of Section IIII-Dll on the typical descending tree. Hence, there 
is hope to progress on this and on related cases. This will 
however not be pursued here and is left for future research. 

V. Numerical Results 

In this section, we study the proposed adaptive spatial 
Aloha schemes quantitatively. We compute various metrics 
formulated through the stochastic geometry based analysis, 
and we also perform simulation. The simulation not only 
validates the analytical model, but also illustrates the perfor- 
mance of schemes for which we do not have an analytical 
characterization. 

A. Computation of the Integrals 

We used Maple and Matlab to evaluate the integrals of 
Section |IVl The infinite integral that shows up in the ex- 
pression of the Laplace transform (flU is handled without 



truncation by Maple and Matlab. The singularity at «; = 
in the contour integrals (fTZt leveraging Parseval's theorem is 
a false singularity and it is also handled without further work 
by either Maple or Matlab. The Matlab code is particularly 
efficient and is used throughout the analytical evaluations 
described below. 

B. Simulation Setting 

We consider a two dimensional square plane with side 
length L, and A'^ nodes placed independently over the plane 
according to the uniform distribution; this corresponds to 
A = N/L^ in the stochastic geometry model0 Each node has 
its receiver randomly located on the unit circle around it, again 
as per the uniform distribution. Thus ra — 1 for all i. We 
set Q = 4 and T = 10. To nullify the edge effect, we take 
into account only the nodes falling in the L/2 x L/2 square 
around the center while computing various metrics. While 
all other parameters remain, we vary L and for different 
simulations. For each parameter set we calculate the average 
of the performance metric of interest over 1000 independent 
network realizations. 

C. Joint Validation of the Analysis and the Simulation 

We validate the analytical expression against the simulation 
for the case of proportional fair medium access. For illustra- 
tion, we plot the cumulative distribution function (c.d.f.) of the 
MAP in Figure ID Here we set L = 40 and consider two values 
of iV, A'' = 400 and = 800, which correspond to A = 0.5 
and A = 0.25 respectively. The plots show that the stochastic 
geometry based formula (see Theorem 14. 2) quite accurately 
predicts the nodes' behavior in simulation. 




simulation 
analysis 
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c.d.f. of MAP 

Fig. 2. Cumulative distribution function of tlie MAP for tlie proportional 
fair case. 

We also study the distributions P°'^(p2; > p) (also referred 
to as fr for 1 1 a; 1 1 = r) defined in Theorem 14.3 1 Figure |3] shows 

^A finite snapshot of a Poisson random proces would contain a Poisson 
distributed number of nodes. However, for lai'ge AL^, the Poisson random 
variable with mean is highly concentrated around its mean. Thus we 
can use nodes for all the realizations in our simulation. 



their plots for A = 0.25 and two values of ||x||, ||a;|| = 1 a 
||x|| = 10. Again the plots based on the analytical expressi 
and those based on simulation closely match. Clearly, unc 
P", nodes closer to the origin are more likely to be inactiw 
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Fig. 3. Demonstration of f^, the MAP distribution of a node at distance r 
from origin, under P". 



D. Performance of the Adaptive Protocols 

In this section we illustrate the performance of various 
adaptive schemes and their benefit over plain Aloha. We com- 
pute the performance metrics via simulation and also through 
analytical expressions whenever the latter are available. In 
such cases the analytical results and the simulation validate 
each other. 

First we set L = 20 and A*' = 50. We consider the maximum 
throughput medium access, however, only accounting for the 
closest interferers. In figure 2] we show steady state behavior 
of our Gibbs sampling based algorithms; we have set the 
temperature r(t) = l/log(l + t). As expected, the improved 
maximum throughput medium access (see the discussion at 
the end of Section IIII-Dll ) insures a better exclusion behavior 
Under this scheme, a lesser number of nodes transmit, and 
neighboring nodes are unlikely to transmit simultaneously. So, 
this is expected to deliver better aggregate throughput. 

Now we keep L fixed at 20, but vary from 10 to 
100; this corresponds to varying A from 0.025 to 0.25 in the 
analytical expressions. We evaluate the aggregate throughputs 
of various Aloha schemes including plain Aloha. The average 
throughputs are plotted in Figure |5] Although in some of 
the schemes we derive the attempt probabilities only con- 
sidering the closest interferers, we always take into account 
the aggregate interference while calculating the throughput. 
When the number of nodes is small, both the throughput 
maximizing medium access and plain aloha have identical 
performance; both prescribe attempt probabilities close to one 
for all the nodes. When the number of nodes increases beyond 
45, the throughput maximizing medium access significantly 
underestimates the interference, and thus its performance de- 
teriorates. On the other hand, the aggregate interference based 
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Fig. 4. Ttirougliput maximizing medium access: Tliere are 100 transmitter re- 
ceiver pairs. The diamonds represent candidate transmitters, and the connected 
circles the corresponding receivers. The solid diamonds represent the nodes 
that transmit in all the slots; others never transmit. The left plot corresponds 
to the maximum throughput medium access and the right one to its improved 
version. 



proportional fair scheme significantly outperforms plain Aloha 
in terms of aggregate throughput also. This benefit is sustained 
even as the number of nodes increases. We also notice that the 
unproved version of maximum throughput medium access (see 
the discussion at the end of Section IIII-Dll ) yields best 
performance among all the schemes, and its performance does 
not deteriorate until a much higher number of nodes. 
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Fig. 5. Thi'oughputs of various medium access schemes as a function of 
the number of nodes. MT, PF, CI and AI stand for maximum throughput, 
proportional fair, closest interferer and aggregate interference respectively. 

In Figure |6] we plot the logarithms of the aggregate 
throughputs corresponding to the two proportional fair medium 
access schemes. The figure illustrates that as the number 
of nodes increases, the performance of the closest interferer 
based medium access worsens in comparison to the perfor- 
mance of the aggregate interference based scheme - this is 
not visible merely looking at the corresponding aggregate 
throughputs (see Figure |5]l). 

In Figure |7] we plot the c.d.f. of the MAP for the propor- 
tional fair case. Our objective is to compare the case when 
nodes account for the aggregate interference with when they 
account for the closest interferer only. We set L — AO and 
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Fig. 6. Performance of the two proportional fair medium access sciiemes as 
a function of the number of nodes. 

plot MAP distributions corresponding to two values of A^, 400 
and 800 (corresponding to A = 0.25 and A = 0.5 respectively). 
As expected, the nodes attempt more aggressively when they 
account for the closest interferer only. While there is a favored 
probability interval in the aggregate interference case, there 
are more than one such intervals in the closest interferer case. 
Also, in the latter case, about 34% of nodes attempt in almost 
all slots, irrespective of N, the total number of nodes. This 
can be understood by noticing that the probability that a node 
is not the closest interferer to any other node is not sensitive 
to N. 
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Fig. 7. Cumulative distribution function of the MAP for the proportional 
fair case. 



VI. Conclusion 

We have shown the feasibility of the performance analysis 
of distributed adaptive protocols that aim at maximizing some 
global utility in a large random network using stochastic 
geometry. 



More precisely, the most natural distributed adaptation of 
the medium access probability of Aloha that aims at pro- 
portional fairness optimization was shown to have a tractable 
optimal MAP distribution. This distribution is obtained from 
the law of a certain shot noise field that describes the in- 
terference created by a typical node to all receivers but his. 
In the Poisson case, the distribution of the optimal MAP is 
obtained as a non-singular contour integral which is amenable 
to an efficient evaluation using classical numerical tools. The 
network performance at optimum can in turn be deduced from 
the latter using Campbell's formula. 

This approach is shown to provide an analytic way of quan- 
tifying the gains brought by this proportionally fair adaptive 
version of Aloha compared to plain Aloha. 

This line of thoughts opens several research directions. The 
first one is the extension to other types of fairness, still in 
the framework of Aloha. We indicate that this is possible at 
least under certain simplifications of the interference model. 
The second and broader question is whether this approach 
can be extended to MAC protocols other than Aloha. An 
example would be an adaptation of the exclusion radius of 
CSMA/CA to the full environment of a node aiming at 
maximizing some utility of the throughput. A third general 
question concerns the evaluation of the "price of decentraliza- 
tion". When the discussed protocols are suboptimal because of 
their greedy/distributed nature, is it possible to use stochastic 
geometry to evaluate the typical discrepancy between the per- 
formance of the distributed scheme and the optimal centralized 
one? 
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Appendix 

Let * = {Y{n} be a point process with marks and V ° be 
the point process V shifted by —z. Given that Yq = {ro,4>), 

p„ = h{\\Xn - {ro, 0)11, \ {{ro,4>),Yn} o ex J, 

where the mapping /i(w, associates with * and the real 
number u the solution of 

1 _ 1 ^ 1 

with the usual convention if there is no solution in [0, 1]. It 
follows from Slivnyak's theorem that 




Pn 



^ '-J 



hi\\Xn-Yo\\,^'-\{Yo,Yn}oex^ 



— T¥S + '- 



1 - H\\Xn-{ro,m,^^ \{Yn}o0xJ 



\\Xn-(ro,. 



+ 1 



It now follows from Campbell's formula that 

^ ^,^°M' ll^n-(ro,^)|h I 1 



log 1 



-I I 



Ja:- (ro,0)|l° +ro 



Ada; 



F^m^-{ro,m,'f'''\{Yo}) = dv) 



vro 



— )Adx 



,|a;- (ro,0)l|° +ro. 



where the last relation follows from Slivnyak's theorem. 



